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L^-REDUCIBILITY AND LOCALIZATION FOR QUASIPERIODIC 

OPERATORS 

SVETLANA JITOMIRSKAYA AND ILYA KACHKOVSKIY 


Abstract. We give a simple argument that if a quasiperiodic multi-frequency Schrodinger 
cocycle is reducible to a constant rotation for almost all energies with respect to the density 
of states measure, then the spectrum of the dual operator is purely point for Lebesgue almost 
all values of the ergodic parameter 9 . The result holds in the setting provided, in addition, 
that the conjugation preserves the fibered rotation number. Corollaries include localization for 
(long-range) ID analytic potentials with dual ac spectrum and Diophantine frequency as well 
as a new result on multidimensional localization. 


1. Introduction 

The relation between spectral properties of discrete quasiperiodic Schrodinger operators and 
their duals, known as Aubry duality, has a long history starting from p]. For nice enough n, 
sufficiently regular normalizable solutions of H'lp = Eip with 

(1.1) {H{x)'lp)n = -Ipn+l + V'n-l + v{x + na)'lpn 
yield Bloch wave solutions of Hip = Eip for 

(1.2) {Smu = E «•” /Ipm-m' + 2 COS 2Tl{am + 9)lpm- 

m'eZ 

where 

v{x) = 

fcez 

Conversely, nice enough Bloch waves yield normalizable eigenfunctions. There was an expec¬ 
tation stemming from [2] that this should translate into the duality between pure point and 
absolutely continuous spectra for H and H. The fact that such correspondence in the direction 
“absolutely continuous spectrum for H implies pure point spectrum for H” could be false even 
for the nicest of v was hrst understood in [27] which was a surprise at the time. Y. Last showed 
that for the almost Mathieu operator Hx dehned by fll.ip with v{x) = 2Acos27ra: and Liouville 
a there is absolutely continuous component in the spectrum of H for A < 1, while H has purely 
singular continuous spectrum. By now it is known that for A < 1 the spectrum of H is purely 
absolutely continuous for all x, a, while for the corresponding H there is an interesting phase 
diagram of singular continuous and pure point phases depending on the arithmetic properties 
of both a and x PIM]. It should be noted that the almost Mathieu family is self-dual 

with Hx = \Hijx- 

A more rehned notion than absolutely continuous spectrum can be given in terms of re- 
ducibility of the corresponding Schrodinger cocycle, see ini. Reducibility, for instance, does 
imply that the generalized eigenfunctions behave as Bloch waves. 

1 




2 


S. JITOMIRSKAYA AND I. KACHKOVSKIY 


The duality has been explored at many levels, from physics motivated gauge invariance back¬ 
ground [29], to operator-theoretic [18], to quantitative [8]. The main recent thrust, however, 
has been in establishing absolute continuity of the spectrum of H [T2], reducibility [30], or more 
subtle property called almost reducibility [H], by proving Anderson localization (or almost local¬ 
ization) for the dual operator H. This was motivated by the development of non-perturbative 
localization methods, starting with [23] , that led to duality-based reducibility conclusions with¬ 
out the use of KAM and thus allowing for much larger ranges of parameters. 

However, recent celebrated results such as Bia provide methods to establish non-perturbative 
reducibility directly and independently of localization for the dual model. Thus, a natural ques¬ 
tion arises: does reducibility of the Schrodinger cocycles imply localization for H{d)l Since nice 
enough reducibility directly yields nicely decaying eigenfunctions for the dual model, the main 
question is that of completeness. This has been a stumbling block in corresponding arguments, 
e.g. [32]. The question is non-trivial even for the almost Mathieu family, where the conjectured 
regime of localization for a.e. x has been < A where jS is the upper rate of exponential 
growth of denominators of continued fractions approximants to a [22], yet direct localization 
arguments have stumbled upon technical difficulties for A < see [SHIE]. At the same 

time, in the recent preprint [9], the Almost Reducibility Conjecture, recently resolved by Avila 
(see [3] where the /3 > 0 case that is needed for [9] is presented), is combined with the tech¬ 
nique from [511201132] to establish reducibility for the dual of the entire < A region, thus 
leading to the above question of completeness of the dual solutions. The authors of [9] use 
certain quantitative information on the model, in particular, the existence of a large collection 
of eigenfunctions with prescribed rate of decay to obtain their completeness, which, ultimately, 
implies arithmetic phase transitions for the almost Mathieu operator. 

In the present paper, we obtain an elementary proof of complete localization for the dual 
model under the assumption of L^-reducibility of the Schrodinger cocycle for H(x) for almost all 
energies with respect to the density of states measure (see Theorem 13.ip . We do not make any 
further assumptions on H{x) or it’s spectra, however, we require the conjugation to preserve 
the rotation number of the Schrodinger cocycle (which is always possible if the conjugation is 
continuous). 

Our result implies a simple proof of localization as a corollary of dual reducibility for the 
almost Mathieu operator throughout the entire < A region (first proved in [9]). Also, if 
the spectrum of H{x) is purely absolutely continuous, v is analytic, and a is Diophantine, the 
assumptions of Theorem 13.11 automatically hold due to [S], which gives a new interpretation 
of duality in the direction from absolutely continuous spectrum to localization, as originally 
envisioned. The argument applies to the multi-frequency case, thus allowing some conclu¬ 
sions on multidimensional localization (see Theorem 13.4p . It has further implications to the 
quasiperiodic XY spin chain, to be explored in a forthcoming paper [2T] . 

2. Preliminaries: Schrodinger cocycles, rotation number and duality 
Let 

{H{x)i))n = i^n+l + + v{x + na)tpn, 

where n G Z, x = (xi,..., x^), a = (ai,..., ad) G T'’*, and a is incommensurate in the sense that 
{1, «!,..., ad} are linearly independent over Q. The family {H{x)}x^fd is called a d-frequency 
one-dimensional quasiperiodic operator family. We will identify T'’* with [0,1)'^, and continuous 
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functions on with continuous Z'^-periodic functions on The eigenvalue equation 
(2.1) ijn+l + i^n-l + v{x + na)'tpn = 

can be written in the following form involving transfer matrices, 


'^n 

i’n-l 




where 


Sv,e{x) = 


E 


vix 


0 


and the pair (a, S^^e) is a Schrodinger cocycle understood as a map (a, Sy^E) : T'^xC^ —)■ T'^xC^ 
given by {a,Sy^E) ■ {x,w) i—)■ (x + a,Sy^E{x) ■ w). Replacing Sy^E with A e SL(2,R) gives a 
dehnition of an SL(2, R)-cocycle. 

Suppose that A is an SL(2, R)-cocycle homotopic to the identity (for example, a Schrodinger 
cocycle). Then there exist continuous functions tc: x T —)• R and u: x T —>■ R"*" such that 

., f cos27r-u\ . . { cos2'Kiy + w{x,y))\ 

Aix) . o = u(x,y) . o / , / \\ ■ 

^ ' ysm 27 r|/ j '^1^1 ysm 2 TT{y + w{x,y)) J 

Let /i be a probability measure invariant under transformation of x T defined by [x, y) t-)- 
{x + a,y + w{x,y)) such that it projects to Lebesgue measure over the first coordinate. Then 
the fibered rotation number of {a, A) is defined as 


p{a,A) = 


'jd-xT 


w dp ) mod Z. 


\i A = Sy^Ei then we will denote its rotation number by p{E), ignoring the dependence on a 
and V. 

If H{-) is a quasiperiodic operator family, and Ee(x){^) are the corresponding spectral pro¬ 
jections, then the density of states measure is defined as 

N{A)= [ {EHix){A)6o,6o)dx, 

JTd 

A C R is a Borel set. The distribution function of this measure N{E) = N{{—oo,E)) = 
N{{—oo,E]) is called the integrated density of states (IDS). A remarkable relation (see, for 
example, [25l [10] for one-frequency case and [16] for the general case) is that for Schrodinger 
cocycles we have N{E) = 1 — 2p{E), and p is a continuous function mapping S onto [0,1/2], 
where S = a{H{x)) is the spectrum which is known to be independent of x. The following 
holds 

Proposition 2.1. Suppose that E ^ T,. Then 2p(E) G OiZ -f 02^ -|- ... -f o^Z -|- Z. 


Aubry duality is a relation between spectral properties of H{x) and the dual Hamiltonian 
(Ol). We will formulate it in the multi-frequency form, with the dual operator fll.2p now being 
an operator in 

I'ljjm-m' + 2 COS 27r(a • m -|- 9)lprm 
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where 

v{x) = 

fcez 

Denote the corresponding direct integral spaces (for H and H respectively) by 


/*© /*© 

s) = P{z)dx, f{z^)de. 

Jj 

Consider the unitary operator U: S) ^ Sj defined on vector functions d' = d' (a:, n) as 

m) = ^(m, 6 -\- a ■ m), 

where denotes the Fourier transform over x G —)■ m G combined with the inverse 

Fourier transform n G Z —)■ 6^ G T. Let also 

/*© /*© 

n:= H{x)dx, ■H:= H{e)de. 

Jfd JTT 

Aubry duality can be formulated as the following equality of direct integrals. 

( 2 . 2 ) UHU-^ = H. 


It is well known (see, for example, [18]) that the spectra of H{x) and H{9) coincide for all x, 9. 
We denote them both by E. Moreover, the IDS of the families H and H also coincideEl 

A continuous cocycle {a, A) is called L'^-reducible if there exists a matrix function B G 
L^(T'^; SL(2, M)) and a (constant) matrix A* G SL(2,M) such that 

(2.3) B{x + a)A{x)B{x)~^ = A*, for a. e. x G T'^. 


We will call a cocycle A -degree 0 reducible, if fl2.3p holds with (possibly complex) B with 
I det B{x)\ = 1 and 


(2.4) 


/g2^*p(E) Q \ 

0 ) ■ 


3. Obtaining localization by duality 
Our main result is as follows. 

Theorem 3.1. Suppose that H{x) is a continuous quasiperiodic operator family such that for 
almost all E E T, with respect to the density of states measure the cocycle is B^-degree 0 
reducible. Then the spectra of dual Hamiltonians H{9) are purely point for almost all 9 eT. 

The following fact is standard and shows that the notion of degree 0 reducibility is needed 
only in discontinuous setting. 

Proposition 3.2. Suppose {a, A) is a continuously reducible cocycle such that fl2.3p holds with 
A* being a rotation matrix. Then it is B^ (and continuously) degree 0 reducible. 

^While the proof in |18] is claimed for a more specific one-frequency family, it applies without any changes 
in the mentioned generality. 
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Proof. Assume that fl2.3p holds with B continuous and homotopic to a constant matrix. Since 
such conjugations B preserve rotation numbers, the matrix A* will be a rotation by the rotation 
number of {a, A), and fl2.4|) can be obtained by diagonalizing A* over C (the diagonalization is 
another conjugation with a constant complex matrix B). 

In general, any continuous map from i?; —)■ SL(2,M) is homotopic to a rotation R^.x for 
some n G lA. If i?; —)■ SL(2,M) satisfies fl2.3p and is homotopic to Rnx, then R_n.xB{x) is 
homotopic to a constant and also satisfies fl2.3p . ■ 

Corollary 3.3. Suppose that v is real analytic and d = 1, a is Diophantine, and that the 
spectrum of the original operator fll.ip is purely absolutely continuous for almost all x. Then 
the spectrum of the dual operator (II■2p is purely point for almost all 6, and the eigenfunctions 
decay exponentially. 

Proof. From |5l Theorem 1.2], it follows that for Lebesgue almost all E E the cocycle (a, 
is rotations reducible. If a is Diophantine, it will also be (analytically) reducible to a constant 
rotation (see, for example, [30]). Hence, we fall under the assumptions of Theorem 13.11 ■ 

We can also formulate the following (and now elementary) result on multi-dimensional local¬ 
ization. 

Theorem 3.4. If a is Diophantine, then there exists Ao(a,n) such that the long-range lattice 
operator 

{H{9)'ljj)m= ^ Vm'Am-m' +‘2X cos 271 {a ■ m + 9)^prn 
has purely point spectrum for X > Ao(q!, v) and almost all 6^ G T. 

Proof. By Eliasson’s perturbative reducibility theorem (see mm]), the cocycle (a, Sx-iy^E) is 
continuously (even analytically) reducible for A > Ao(a,n) for almost all values of the rotation 
number. Hence, it satisfies the hypothesis of Theorem 13.11 ■ 

Remark 3.5. We are not aware of other multidimensional localization results with purely 
arithmetic conditions on the frequency. In [13], Theorem 13.41 was proved for general poten¬ 
tials with two critical points (rather than for cosx). Long range operators were also considered. 
However, the result was obtained for an unspecified full measure set of frequencies. In the book 
P, a possibility of alternative proof was mentioned (for general analytic potentials), with a 
part of the proof only requiring the Diophantine condition on a. 

Remark 3.6. For every v and for almost every E among those for which the Lyapunov exponent 
of the Schrodinger cocycle (a, Sx^e) is zero, this cocycle can be L^-conjugated to a cocycle of 
rotations R 4 >{x), see [26][H], or [3 Theorem 2.1] with a reference to [31]. Further conjugation to 
a constant rotation, such as in the assumption of Theorem 13.11 requires solving a cohomological 
equation (p{x + a) — (p{x) = (j){x) modZ. 

In the remaining two sections we prove Theorem 13.11 

4. Covariant operator families with many eigenvectors 

In this section, we make the following observation: if, for a covariant operator family {H{u)}, 
we have found a single eigenvector of H{u}) for each u, and the eigenvectors for different uj on the 
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same trajectory are essentially different (that is, not obtained from each other by translation), 
then, for almost every ca, the eigenvectors of H (ca) obtained by such translations form a complete 
set. 

Let (n, du) be a measurable space with a Borel probability measure, equipped with a family 
{T^,n G of measure-preserving one-to-one transformations of hi such that 
Let {H{u;),u; G hi} be a weakly measurable operator family of bounded operators in 
with the property iL(T”a;) = TnH{uj)T-n, where Tn is the translation in Z'^ by the vector n. 
We call {H{u)} a covariant operator family with many eigenvectors if 

(1) there exists a Borel measurable function i?: hi —)■ M such that for almost every u E Q 
the equality E(T^u)) = E(T^u;) holds only for m = n. B 

(2) There exists a function m: hi —)■ E {!/'), not necessarily measurable, such that for almost 
every cj G hi we have H{u)u{u) = E{u)u{u) and ||M(a;)||p(zd) = 1 . 

Theorem 4.1. Suppose that {H{u)} is a covariant operator family in Z'^ with many eigenvec¬ 
tors. Then, for almost every u, the operator H{u) has purely point spectrum with eigenvalues 
E(T~‘^u) and eigenvectors Tnu{T~‘^u). 

Proof. Let Ekipj) := E{T~^u), and Pkipj) be the spectral projection of H{uj) onto the eigenspace 
corresponding to Ek{uj). Both these functions are weakly measurable in u. For almost all u, 
we have for all k 

H{T^ijj)Tku{oj) = TkH{uj)u{uj) = E{oj)Tku{oj) = Ek{T^uj)Tku{uj). 

Hence, Tkuioj) belongs to the range of PkiT^oj), and, for each basis vector 5i G we have 

( 4 . 1 ) {Pk{T^u)5u5i) ^ \{Tku{u),5i)\\ 

Let P{uj) := By assumption, all Ek{uj) are different (for almost all uj), and hence 

P{uj) is a projection. We have 


(4.2) [{P{oj)6i,6i)doj = f '^{Pk{uj)Si,Si) duj = 

Q i 

= f Y.^Pk{T^uj)5u5i)du^ I J2\iTku{u),6i)\‘^du = l. 

i i keZ'i 

Since the left hand side is bounded by 1, the last inequality must be an equality, and hence, 
for any /, the inequality fl4.1l) must also be an equality for almost every cu. By passing to the 
countable intersection of these sets of u, we ultimately obtain a full measure set for which 
P{(.v) = /, and hence the eigenfunctions from the statement of the theorem indeed form a 
complete set. ■ 


Remark 4.2. The left hand side of (14. 2 h is the “partial density of states measure” of hi in the 
sense that it is the density of states obtained only from a part P{oj) of the spectral measure. The 
fact that this measure is still equal to 1 implies that, for almost every u, that part was actually 
the full spectral measure, which is equivalent to completeness of the eigenfunctions. The result 
also implies simplicity of the point spectrum and that the function u{-) can actually be replaced 

^In other words, E should be one-to-one on almost all trajectories. However, the wording “almost all 
trajectories” is ambiguous as there is no natural measure on the set of trajectories. 
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by a measurable function with the same properties (because spectral projections are measurable 
functions, and, due to simplicity, one can choose measurable branches of eigenfunctions; see 
[T5] for details). 


5. Proof of Theorem 13.11 

Let ©jf be the set of all 9 G [0,1/2] such that p{E) = 9 for some E E T, with being 
L^-degree 0 reducible, and 29 ^ aiZ + ... + + Z. By the assumptions, this set has full 

Lebesgue measure in [0,1/2]. Moreover, for 9 G ©j!" there exists exactly one E such that 
p{E) = 9, because p is nonincreasing and can be constant only on intervals contained in M \ S, 
on which the rotation number must be rationally dependent with a, see Proposition 12.11 Let 
us denote this inverse function by E{9), initially defined on ©+, extend it evenly onto —©jf, 
and then 1-periodically onto M. Let us now consider 

fir 0 ) = /Bri'(i^..E(«))/i|Bri‘(-.B(»))iu»(T). 0 6 |o,i/ 2 ] 

’ \B:px,E(-e))/\\Br,'(:,E(-e))h,(r,, 9 e 1-1/2,0), 

and then also extend f{x, 9) 1-periodically over 9. Denote the domain in 9 of both functions 
by ©r which is a similar (even, and then periodic) extension of ©/“. For 9 E Or, fhe function 
/ generates Bloch solutions 

'ilj{x,9)n = f{x + na,9), n G Z, 

which, for almost all x, are formal solutions of the eigenvalue equation {E[{x)'ip)n = E{9)'4)n. It 
is also easy to see that, if u{9)m = /(tu, 9) is the Fourier transform of / in the variable x, then 

{H{9)u{9))^ = E{9)u{9),n, m E Z'". 

Finally, E{9 + k- a) = E{9 + l ■ a) with k ^ I can only happen if 26* G Z -|- Zai Za^^, and 

so, the function E{-) satisfies the assumptions of Theorem 14.11 from which the result follows. ■ 

Remark 5.1. For almost all 9, we, in addition, can choose f{x,—9) = f{x,9). This is a 
standard argument (see, e.g., [30], [0] , or |3] in the discontinuous setting). Let 


E{x,9) 


( fix, 9) fix, 9) \ 

\e“^"®/(x — a, 9) e^’^*®/(x — a,9)J ' 


From the equation, it is easy to see that S^^EEix, 9) = F(x -t- a, 9)A^, and hence (i(x, 9) = 
det E{x, 9) is a function of x invariant under x ^ x + a. Therefore, it is almost surely a 
constant function of x. If dix,9) ^ 0, then we can simply choose B{x) = E[x,9)~^. Suppose 
that (i(x, 9) = 0. Then we have 

fix, 9) = (fix, 9)f{x, 9), e^'^'-^fix + a,9) = p{x, 9)e~‘^^^^fix + a,9). 

for some function (p with |(p(x, 6*)| = 1 and (p(x -f- a, 6*) = e“'^”®(p(x, 9). That is only possible if 
29 is a rational multiple of a modulo Z, which excludes at most countable set of 9s. 
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